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The impact of the experimental acceptance, i.e. transverse-momentum (pr) cut-off and limited 
rapidity region, on the earlier predicted irregularity in the excitation function of the baryon stop¬ 
ping is studied. This irregularity is a consequence of the onset of deconfinement occurring in the 
compression stage of a nuclear collision and manifests itself as a wiggle in the excitation function 
of the reduced curvature (C y ) of the net-proton rapidity distribution at midrapidity. It is demon¬ 
strated that the wiggle is a very robust signal of a first-order phase transition that survives even 
under conditions of a very limited acceptance. At the same time the C y for pure hadronic and 
crossover transition scenarios become hardly distinguishable, if the acceptance cuts off too much 
of the 1 ow-pt proton spectrum and/or puts too narrow rapidity window around midrapidity. It is 
found that the shape of the net-proton rapidity distribution near midrapidity depends on the pr 
cut-off. This implies that the measurements should be taken at the same acceptance for all collision 
energies in order to reliably conclude on the presence or absence of the irregularity. 
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I. INTRODUCTION 

The onset of deconfinement in relativistic heavy-ion 
collisions is now in the focus of theoretical and exper¬ 
imental studies of the equation of state (EoS) and the 
phase diagram of strongly interacting matter. This prob¬ 
lem is one of the main motivations for the currently 
running beam-energy scan [lj] at the Relativistic Heavy- 
Ion Collider (RHIC) at Brookhaven National Laboratory 
(BNL) and the low-energy-scan program @ at Super Pro¬ 
ton Synchrotron (SPS) of the European Organization 
for Nuclear Research (CERN) as well as for construct¬ 
ing the Facility for Antiproton and Ion Research (FAIR) 
in Darmstadt [sj and the Nuclotron-based Ion Collider 
Facility (NICA) in Dubna Q. 

In Refs. OH it was argued that the baryon stopping 
in nuclear collision can be a sensitive probe for the onset 
of deconfinement. Rapidity distributions of net-protons 
were calculated 0, §] within a model of the three-fluid 
dynamics (3FD) [9f in scenarios with and without decon¬ 
finement transition. These calculations were performed 
employing three different types of EoS: a purely hadronic 
EoS pj| (hadr. EoS) and two versions of the EoS involv¬ 
ing deconfinement 11]. The latter two versions are an 
EoS with a first-order phase transition (2-phase EoS) and 
one with a smooth crossover transition (crossover EoS). 

It was found that 3FD predictions within the first- 
order-transition scenario exhibit a “peak-dip-peak-dip” 
irregularity in the incident energy dependence of the form 
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of the net-proton rapidity distributions in central colli¬ 
sions. At low energies, rapidity distributions have a peak 
at the midrapidity. With the incident energy rise it trans¬ 
forms into a dip, then again into a peak, and with further 
rising energy the midrapidity peak again changes into a 
dip, which already survives up to arbitrary high ener¬ 
gies. The behaviour of the type “peak-dip-peak-dip” in 
central collisions within the 2-phase-EoS scenario is very 
robust with respect to variations of the model parame¬ 
ters in a wide range. This behaviour is in contrast with 
that for the hadronic-EoS scenario, where the form of 
distribution at midrapidity gradually evolves from one 
with a peak to one with a dip. The case of the crossover 
EoS is intermediate. Only a weak wiggle of the type of 
“peak-dip-peak-dip” takes place. 

Experimental data also reveal a trend of the “peak- 
dip-peak-dip” irregularity in the energy range 8A GeV 
< Ei a i, < 40 A GeV, which is qualitatively similar to that 
in the first-order-transition scenario while quantitatively 
it differs. It is very likely that the quantitative discrep¬ 
ancy is due to the inadequacy of the model for the 2-phase 
EoS for which the onset of the phase transition lies at 
rather high densities, above 8 times saturation density 
at low temperatures. Recent models suggest an onset 
density of half that value |~i~2l IT~il which would place the 
onset of the wiggle structure closer to the experimentally 
indicated position. However, the experimental trend is 
based on preliminary data at energies of 20A GeV and 
30A GeV. Therefore, updated experimental results at en¬ 
ergies 20A and 30A GeV are badly needed to pin down 
the preferable EoS and to check the hint to the wiggle be¬ 
havior of the type “peak-dip-peak-dip” in the net-proton 
rapidity distributions. Moreover, it would be highly de¬ 
sirable if data in this energy range were taken within the 
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same experimental setup and with the same experimental 
acceptance. 

The calculations of Refs. 00 were performed assum¬ 
ing the acceptance for net-protons to be wide enough to 
include almost all emitted particles. In practice, the ex¬ 
tension of the acceptance beyond the range of 0 < pr < 2 
GeV/c does not practically change the net-proton rapid¬ 
ity distributions in the incident-energy range of interest. 
However, the actual experimental acceptance can be nar¬ 
rower. For the NICA MPD experiment it is restricted by 
the proton identification capabilities in the TOF detector 
to 0.4 GeV/c < pr < 1.0 GeV/c in the central rapidity 
range \y\ < 0.5 E3- For the case of the STAR beam 
energy scan the analysis of the excitation function of the 
net proton rapidity distribution is still under way. In 
this case the acceptance will be restricted to the range 
0.4 GeV/c < p T < 3.0 GeV/c and \y\ <0.5 0. 

The purpose of the present paper is to investigate the 
question whether the “peak-dip-peak-dip” irregularity 
survives when experimental circumstances force a nar¬ 
rowing of the accessible acceptance region. 


II. EQUATIONS OF STATE 
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FIG. 1: Pressure scaled by the product of normal nuclear 
density (no = 0.15 fm” 3 ) and nucleon mass (mjv) versus 
baryon density scaled by the normal nuclear density for three 
considered equations of state. Results are presented for three 
different temperatures T = 10, 100 and 200 MeV (from bot¬ 
tom upwards for corresponding curves). 


Figure [T] illustrates the differences between the three 
considered EoS. The deconfinement transition makes an 
EoS softer at high temperatures and/or densities. The 2- 
phase EoS is based on the Gibbs construction, taking into 
account simultaneous conservation of baryon and strange 
charges. However, the displayed result looks very similar 
to the Maxwell construction, corresponding to the con¬ 
servation of just the baryon charge, with the only differ¬ 
ence that the plateau is slightly tilted, which is practically 
invisible. 

As demonstrated in Refs. 0 0 , the deconfinement 
transition in central Au+Au collisions starts at the top 
AGS energies in both cases. It gets practically completed 
at low SPS energies in the case of the 2-phase EoS. In 
the crossover scenario it lasts till very high incident en¬ 
ergies. We would like to note that the density range for 
the 2-phase EoS in Fig. |T] with an onset of deconfine¬ 
ment above n ~ 8 no is rather high. In recent models 
for the spinodal decomposition accompanying the phase 
transition fl~2l the onset of the phase coexistence is 
at about 4 no. This corresponds to a lower limit for the 
onset of the deconfinement phase transition in cold, sym¬ 
metric nuclear matter which was obtained by a hybrid 
EoS model from constraints on the occurrence of quark 
matter phases in massive neutron stars [l5j , see also [l6[ . 


III. IRREGULARITY AT CONSTRAINED 
ACCEPTANCE 

The calculations at all collision energies were per¬ 
formed for Au+Au (6 = 2 fm) central collisions despite 
the fact that some experimental data were taken for cen¬ 


tral Pb+Pb collisions. This was done in order to avoid 
uncertainties associated with different colliding nuclei. 
However, in fact at the same incident energy the com¬ 
puted results for Pb+Pb collisions at 6 = 2.4 fm are very 
close to those for Au+Au at 6 = 2 fm. The calculations 
were performed for four different acceptance ranges for 
the transverse momentum (pt) and rapidity ( y ) of the 
measured proton: 

i. 0 < pr < 2 GeV/c and a very unrestrictive con¬ 
straint to the rapidity range \y\ < 0.7 ybeam, where 
2/beam is the beam rapidity in the collider mode, 
which is practically equivalent to the full accep¬ 
tance; 

ii. 0.4 < pt < 1 GeV/c and \y\ < 0.5, the expected 
MPD acceptance [l7j; 

iii. 1 < pr < 2 GeV/c and \y\ < 0.5, an acceptance 
range where low-momentum particles witnessing 
collective behaviour are largely eliminated; 

iv. 0.4 < pt < 3 GeV/c and \y\ < 0.5, the range of the 
STAR acceptance jl8j . 

We separately study effects of the pr and y constraints 
in order to reveal their relative importance. 

A direct measure of the baryon stopping is the net- 
baryon (i.e. baryons-minus-antibaryons) rapidity distri¬ 
bution. However, since experimental information on neu¬ 
trons is unavailable, we have to rely on net-proton (i.e. 
proton-minus-antiproton) data. Presently there exist ex¬ 
perimental data on proton (or net-proton) rapidity spec¬ 
tra at AGS [311 and SPS [23l - [27t energies. At AGS 






3 



FIG. 2: Midrapidity value of the net-proton rapidity distribution as a function of the collision energy in central (b — 2 fm) 
Au+Au collisions in different windows of the transverse momentum pr'- (a) 0 < pr < 2 GeV/c, (b) 0.4 < pr < 1 GeV/c, (c) 
1 < Pt < 2 GeV/c and d) 0.4 < pr < 3 GeV/c. Results for the wide py-window (0 < pr < 2 GeV/c) are also presented 
in panels (b) - (d) by the corresponding thin lines for the sake of comparison. Results of simulations with different EoS’s are 
presented. Experimental data are from the collaborations E895 [lf|, E877 j2cj, E917 [21j], E866 [22], NA49 [2p| - l27ll and STAR 


energies, the yield of antiprotons is negligible, therefore 
the proton rapidity spectra serve as a good probe of the 
baryon stopping. 

In order to quantify the above-discussed “peak-dip- 
peak-dip” irregularity, it is useful to make use of the 
method proposed in Ref. [H|. For this purpose the data 
on the net-proton rapidity distributions are fitted by a 
simple formula 

cLN 

— = a (exp {—(l/w s ) cosh(y — y s )} 
ay 

+ exp{-(l/w s )cosh(y + y s )}) , (1) 

where a, y s and w s are parameters of the fit. The form 
m is a sum of two thermal sources shifted by ±y s from 
the midrapidity which is put to be y m id = 0 as it is in 
the collider mode. The width w s of the sources can be 
interpreted as w s = (temperature)/(transverse mass), if 
we assume that collective velocities in the sources have 
no spread with respect to the source rapidities ±y s . The 
parameters of the two sources are identical (up to the 
sign of y s ) because only collisions of identical nuclei are 
considered. 

The above fit has been prepared using the least-squares 
method and applied to both, available data and results of 
calculations. The fit was performed in the rapidity range 
\y\ < 0.7 !/beam, where ybeam is the beam rapidity in the 
collider mode. The choice of this range is dictated by 
the data. As a rule, the data are available in this rapid¬ 
ity range, sometimes the data range is even more narrow 
(80A GeV and new data at 158A GeV [27| )■ We apply 
the above restriction in order to deal with different data 
in approximately the same rapidity range. Another rea¬ 
son for this cut is that the rapidity range should not be 


too wide in order to exclude contributions of cold spec¬ 
tators. The fit in the rapidity range |y| < 0.5 j/beam has 
been also done in order to estimate uncertainty of the fit 
parameters associated with the choice of fit range. An 
additional set of fits has been done under the constraint 
\y\ < 0.5. This fit was applied only to the 3FD simula¬ 
tion results. Due to low number of experimental points 
within the |j/| < 0.5 range and their insufficient accuracy 
such a fit gives too large error bars for the deduced pa¬ 
rameters of the fit © at yj snn > 5 GeV, whereas it is 
practically identical to the \y\ < 0.7 j/beam fit at tJsnn < 
5 GeV. 

A useful quantity, which characterizes the shape of the 
rapidity distribution, is a reduced curvature of the spec¬ 
trum at midrapidity defined as follows 



= (ybe a m/w s ) 2 (sinh 2 y s - w s coshy s ) . (2) 

The factor 1/ (ybe&mdN/dy) 0 is introduced in order to 
cancel out the overall normalization of the spectrum. The 
second part of Eq. © presents this curvature in terms 
of the parameters of the fit ©■ The reduced curvature, 
C y , and the midrapidity value, ( dN/dy ) y=0 , are two in¬ 
dependent quantities quantifying the the spectrum in the 
midrapidity range. Excitation functions of these quanti¬ 
ties deduced both from experimental data and from re¬ 
sults of the 3FD calculations with different EoS’s are dis¬ 
played in Figs. [2] and [3] 

In Fig. [2]the midrapidity values of the rapidity spectra 
were taken directly from experimental data and calcu¬ 
lated results. Therefore, only experimental error bars 
are displayed there. The uncertainty associated with the 
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FIG. 3: Midrapidity reduced curvature [see. Eq. <[2j)] of the (net)proton rapidity spectrum as a function of the center-of-mass 
energy of colliding nuclei as deduced from experimental data and predicted by 3FD calculations with different EoS’s: the 
hadronic EoS (hadr. EoS) [10,] (left column of panels), the EoS involving a first-order phase transition (2-ph. EoS, middle 
column of panels) and the EoS with a crossover transition (crossover EoS, right column of panels) into the quark-gluon phase Ha- 
Upper bounds of the shaded areas correspond to fits confined in the region of \y\ < 0.7 j/beam, lower bounds, |y| < 0.5 i/beam- 
Results are presented for four different windows of the transverse momentum pt : 0 < pr < 2 GeV/c (top row of panels), 
0.4 < pr < 1 GeV/c (second row of panels), 1 < pr < 2 GeV/c (third row of panels) and 0.4 < pr < 3 GeV/c (bottom row of 
panels). Results for the py-window of the top panel with experimental data (0 < pt < 2 GeV) are also presented in the lower 
panels with different restricted pr -windows by shaded areas bounded by the corresponding thin lines for the sake of comparison. 
Results of fits within the range of \y\ < 0.5 are displayed by corresponding bold lines. In those cases, when (|j/| < 0.5) fit 
substantially differs from that in the \y\ < 0.9 range, the (|p| < 0.9) results are also displayed by bold double-dash-triple-dotted 
lines. 
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choice of the rapidity range turned out to be the domi¬ 
nant one for the C y quantities deduced from both exper¬ 
imental data and simulation results. Therefore, in Fig. [3] 
results for the curvature C y in the wide rapidity range 
are presented by shaded areas with borders correspond¬ 
ing to the fit ranges \y\ < 0.7 t/beam and \y\ < 0.5 t/beam- 
The C y results in the narrow rapidity range \y\ < 0.5, 
corresponding to the MPD and STAR acceptance, are 
also displayed by bold lines. In order to control the error 
induced by the narrow range \y\ < 0.5, the C y calcu¬ 
lations were also performed in the range \y\ < 0.9. In 
most cases the (|y| < 0.5) and (|y| < 0.9) results turned 
out to be very close to each other. In those few cases, 
when ( \y\ < 0.5) and (|y| < 0.9) results substantially dif¬ 
fer, the ( \y\ < 0.9) results are also displayed in Fig. [3] 
All presently available experimental data shown in the 
leftmost panel of Fig. [2] and in the top panels of Fig. [3] 
approximately correspond to acceptance range (i). 

As seen from Fig. [2j approximately 60% (for differ¬ 
ent collision energies and EoS’s) of the protons produced 
at midrapidity will be covered by the MPD acceptance 
window, while the corresponding coverage of the STAR 
acceptance is ~80%. 

The irregularity in the data is distinctly seen as a 
strong wiggle in the excitation function of C y , see Fig. [3] 
Various data used to deduce C y approximately corre¬ 
spond to the wide acceptance window (i). Of course, this 
is only a hint to an irregularity since this wiggle is formed 
only by preliminary data of the NA49 collaboration. In 
the wide acceptance window (i) (shaded bands in the up¬ 
per raw of panels in Fig. [3]) the C y excitation function in 
the first-order-transition scenario manifests qualitatively 
(though not quantitatively) the same wiggle irregularity 
(middle-column upper-raw panel in Fig. [3]) like that in 
the data fit, while the hadronic scenario produces purely 
monotonous behavior. The crossover EoS represents a 
very smooth transition. Therefore, it is not surprising 
that it produces only a weak wiggle in C y . 

The application of various pp cuts without confining 
the rapidity range (shaded bands in Fig. [3]) does not 
qualitatively change the picture. As seen from Fig. [3] 
the wiggle in the energy dependence of C y is a very ro¬ 
bust signal of the first-order phase transition. It sur¬ 
vives even at very limited pr-acceptance. The important 
difference between different pr-acceptances is that the 
wiggle is completely located in the range of positive cur¬ 
vatures C y (concave shapes of the rapidity distribution 
near midrapidity) at low-p^ acceptance, like the MPD 
one (ii). While for the high-pr acceptance (iii), the wiggle 
entirely lies in the range of negative C y (convex shape). 
Therefore, the name of ”peak-dip-peak-dip” for this ir¬ 
regularity becomes not quite correct. The amplitude of 
the wiggle becomes somewhat weaker when the pt cut¬ 
off from below is too strong, e.g., when 1 < pp < 2 
GeV/c. It is expectable because the wiggle is an effect 
of collective behaviour of the system, in which predomi¬ 
nantly low-momentum particles participate. If these low- 
momentum particles are cut off by the acceptance, the 


collective effects, in particular, the wiggle, become less 
manifested. 

The C v excitation functions for hadronic and crossover 
scenarios are very similar to each other already in the 
wide acceptance window (i) (the upper raw of panels in 
Fig. [3D. The crossover scenario results in a very weak 
wiggle in C y . Under the MPD (ii) and STAR (iv) pp 
cuts the basic features of the C y excitation functions re¬ 
main the same, only the weak crossover wiggle turns out 
to be shifted to the range of higher (the MPD case) or 
lower (the STAR case) C y . At the high-pr cut (iii) the 
qualitative difference between predictions of the hadronic 
and crossover scenarios is practically washed out. This a 
consequence of the above-mentioned lack of collectivity 
in the behavior of the high-momentum particles. 




FIG. 4: Rapidity distribution of net-protons in central (b = 
2 fm) Au+Au collisions at collision energies snn =7.7 GeV 
(a) and 19.6 GeV (b) for the case of the pr range 0 < pr < 
2 GeV/c. 

The actual MPD (ii) and STAR (iv) acceptance condi¬ 
tions also include the restriction |p| < 0.5. Applying this 
restriction in addition to the pp cuts leads to the corre¬ 
sponding bold lines in Fig. [3] In order to study how re- 
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strictive the condition \y\ < 0.5 is, calculations were also 
done for the rapidity window \y\ < 0.9. In most cases the 
C v results for the range \y\ < 0.5 practically coincide with 
those in the \y\ < 0.9 window and are quite close to those 
in the wide rapidity window \y\ < 0.5 j/beam- In these 
cases the \y\ < 0.9 results are not displayed in Fig. 0 In 
a few cases the \y\ < 0.5 and \y\ < 0.9 results consider¬ 
ably differ. These cases correspond to the hadronic EoS 
when \o~w~pt protons are included in the analysis. For 
these cases the \y\ < 0.9 results are displayed by bold 
double-dash-dash-triple-dotted lines in Fig. [3] (contrary 
to solid bold lines for the \y\ < 0.5 results). 

As seen from Fig. 02 under the constraints of the MPD 
(ii) and STAR (iv) acceptance conditions the wiggle in 
the energy dependence of C y is very robust for the first- 
order phase transition and the crossover one. However, 
this is not the case for the hadronic EoS. The hadronic- 
EoS C y excitation functions now exhibit a local weak wig¬ 
gle in the region 8 < tJsnn <12 GeV (in all cases except 
for that with the high-px range 1 < pr < 2 GeV/c) quali¬ 
tatively similar to that for the crossover transition. More¬ 
over, the C y curvature becomes negative at < y/swN = 
19.6 GeV (again in all cases except for that with the high- 
Pt range). These peculiarities are consequences of a fine 
structure of the rapidity distribution near midrapidity 
that becomes dominant in the narrow rapidity window 
\y\ < 0.5. This situation is illustrated in Fig. [I] It is 
seen that at tJsnn = 7.7 and 19.6 GeV there are tiny 
maxima at the midrapidity in the hadronic-EoS rapidity 
distributions, which result in the weak wiggle and neg¬ 
ative curvature, respectively, observed in Fig. [3] The 
global shape of the hadronic-EoS rapidity distributions 
does not exhibit such peculiarities, as it is clear from C y 
calculations in a wider rapidity range. For a comparison, 
see also Fig. 4 of Ref. j29|. 

Therefore, we can conclude that under the conditions 
of the the MPD and STAR acceptance it is possible to 
distinguish the the first-order phase transition, the onset 
of which is signalled by a strong wiggle in the excita¬ 
tion function of C y . However, the difference between the 
purely hadronic case and that with the smooth crossover 
transition becomes ambiguous. 


IV. CONCLUSION 

An irregularity in the baryon stopping is a natural con¬ 
sequence of deconfinement occurring in the compression 
stage of a nuclear collision. It is a combined effect of 
the softest point {DJ [3l| of an EoS and a change in the 
nonequilibrium regime from hadronic to partonic one. As 
was demonstrated in Refs. [7, §j, this irregularity man¬ 
ifests itself as a wiggle in the excitation function of a 
reduced curvature ( C y ) of the net-proton rapidity dis¬ 
tribution at midrapidity. These calculations were per¬ 
formed in the full acceptance range for protons. In the 
present paper we studied the effect of a restricted accep¬ 
tance, in particular, the one expected for the NICA MPD 


experiment and the one of the STAR beam energy scan 
program, on the earlier predicted C y wiggle. 

It was found that the wiggle in the excitation function 
of C y is a quite robust signal of the onset of deconfine¬ 
ment as a first-order phase transition. It survives even 
at very limited acceptance. Therefore, the MPD experi¬ 
ment will be suitable for the experimental investigation 
of this effect. It was also found that the shape of the 
net-proton rapidity distribution near the midrapidity de¬ 
pends on the experimental acceptance. Hence, only if the 
measurements are taken at the same acceptance for all 
collision energies we can reliably conclude on the pres¬ 
ence or absence of the wiggle. Notice that this is not 
the case for the data analysed so far. These data were 
measured in four different experiments. Therefore, the 
question if the presently available data really indicate a 
wiggle in the C y excitation function is still open. It is 
the more so because some of the data in the region of the 
expected wiggle still have a preliminary status. Data in 
preparation by the STAR experiment will meet the re¬ 
quirement that the excitation function of the curvature 
C y shall be measured by one experiment and under same 
acceptance conditions since they stem from a collider ex¬ 
periment with an acceptance range widely independent 
of the beam energy. These data, however, come from the 
energy range above yj snn = 7.7 GeV, while the wiggle is 
expected in the range between 4 and 8 GeV. Therefore, 
systematic measurements with the NICA MPD collider 
experiment in the whole energy range of interest would 
be highly desirable to clarify this problem. 

It is important to emphasize that the discussed irreg¬ 
ularity is a signal from the hot and dense stage of the 
nuclear collision. It is formed at the (nonequilibrium) 
compression stage of the collision, as it was argued in 
Refs. §]. Therefore, a theoretical approach should 
be able to treat differences in the models for the EoS 
(in particular, absence or presence of a first-order decon¬ 
finement transition) already at this early nonequilibrium 
compression stage of the collision. The 3FD model does 
it by means of the three-fluid interactions. This, at least, 
makes the baryon stopping EoS dependent. At the same 
time, a hybrid fluid model does not distinguish different 
EoS models at the early collision stage because the ini¬ 
tial state for the hydrodynamical evolution is prepared by 
means of the same hadronic kinetic model for any EoS. 
Hence, the hybrid fluid model mainly reveals the baryon 
stopping inherent in the hadronic kinetic model and can¬ 
not produce an irregularity of this stopping relevant to 
the hydrodynamical EoS. Here the situation is similar to 
the case of the directed flow that is also predominantly 
formed at the early nonequilibrium stage of the collision 
(see discussion in Ref. [32jh 
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